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Abstract 

The insulator-metal transition volumes for BaS, BaSe, and BaTe are 
calculated for the first time, using the self-consistent augmented- 
spherical-wave technique. The metallization transition volumes are 1) 
smaller than those corresponding to the NaC £-*CsC i structural transitions, 
but, 2) 10-15% larger than those obtained by the Herzfeld dielectric 
theory. The calculated equilibrium energy gaps in the NaCi structure 
underestimate the measured ones by 50-60%. 
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I. Introduction 

The availability of diamond anvil pressure cells and synchrotron X-ray 
radiation permit precise optical and structural studies of the pressure- 
induced metallization of closed-shell compounds. The metallization is 
usually described 1 within the one-electron approximation as the broadening 
and subsequent overlap of a filled valence band and an empty conduction 
band. The closed-shell Ba chalcogenides are particularly promising 
compounds for metallization: 1) the fractional compression for metalli- 

zation should be less than for isoelectronic compounds and elements, such 
as Csl and Xe, which have larger band gaps, and 2) the metallization 
pressure should be less than for other divalent chalcogenides (e.g., MgTe) 
which have larger bulk moduli. 

In this paper we describe self-consistent band structure calculations 
for BaS, BaSe, and BaTe under compression. The calculations predict band 
overlap transition volumes, illustrate the chemical trends in these 
volumes, and will provide a test of the one-electron approximation used to 
calculate the band gaps when experimental data on the metallization volumes 
becomes available. All of the compounds are studied in both the equil- 
ibrium NaCt(Bl) crystal structure and the CsCz(B2) crystal structures, 
because they are known 2 * 3 * 4 to transform under pressure into the CsCi 
structure. This is the first systematic study of the insulator-metal tran- 
sitions in this series of compounds, taking into account the effects of the 
structural transitions. The volume-dependent band gap of BaO has pre- 
viously been calculated 5 by the APW method, only in the NaC £ structure. 

BaO is not included here because our code does not work correctly for the 
highly compressed Ba atom close to the band overlap transition. 
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Furthermore, accurate treatment of strongly localized orbitals such as 0 2p 
is difficult within density functional theory. 

II. Method of Calculation 

The calculations are performed using the self-consistent aug- 
mented-spheri cal -wave method. 6 In this method the one-electron Schrodinger 
equation is solved approximately using an energy-independent basis set 
obtained by numerical integration of the Schrodinger equation inside atom- 
centered spheres, with boundary conditions which take into account the full 
symmetry of the crystal. The self-consistent one-electron potential is 
generated with a local -density approximation 7 for exchange and correlation 
of the Hedin-Lundqvist form. 8 Relativistic effects are neglected. We now 
comment briefly on the accuracy of the two primary approximations made in 
the calculations: 

A) The local -density approximation (LDA). It is well known tht 
for many closed-shell systems the local -density approximation yields theo- 
retical zero-pressure band gaps which are smaller than the experimental 
gaps by up to a factor of two. In systems with bands that are very narrow' 
compared to the band gap it is possible to remove a large part of the dis- 
crepancy by the use of the self-interaction correction (SIC) scheme. 9 
Unfortunately it is not known how to extend the SIC scheme accurately and 
unambiguously to problems in which the bands are broad compared to the band 
gap, such as the metallization process of interest here. 10 On the other 
hand, since the local -density approximation is expected to improve as the 
system becomes more metallic or free-electron like, it will be more 
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accurate close to the'insulator-metal transition than at the equilibrium 
volume. The discontinuity 11 * 12 between the exact potentials for valence 
and conduction band electrons should decrease with the band gap, and in 
fact vanishes 12 in a simple model calculation. Furthermore, recent 
exchange-only calculations 13 for atoms have shown that the differences 
between eigenvalue spacings obtained using the local density potential and 
the exact ground state potential are small. Thus a fairly accurate 
estimate of the insulator-metal transition volume may result even if the 
zero-pressure band gap is inaccurate. 

B) The use of the ASW basis set. We use the following basis 
orbitals: Ba 6s5p5d4f, S 3s3p3d4f, Se 4s4p4d4f, and Te 5s5p5d4f. In our 
ASW calculations 14 for fee Xe, and Csl in the CsCa structure, we obtain 
gaps which differ from APW results 15 * 16 by less than 0.2 eV. The ASW 
method is also similar in accuracy to the linear muffin-tin orbital (LMTO) 
method with the "combined correction" term. 17 The band structure of Si in 
the diamond 18 and fee 19 crystal structures has been calculated with the 
LMTO method. In the diamond structure, empty spheres were included, 
resulting in eightfold coordination of the spheres. Comparison of the LMTO 
results with full -potential linear-augmented-plane-wave method 20 and 
self-cmsistent pseudopotential 1 9 results again reveals discrepancies of 
less than 0.2 eV. Thus the use of the ASW method for solving the 
one-electron . part of the problem does not appear to cause major errors in 
closely packed structures such as CsCa. However, the NaCt structure 
results are expected to be more sensitive to the approximations in the ASW 
method, because of the lower coordination number. The sensitivity tests to 
be presented in Section III confirm these expectations. 
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III. Results 

The relevant particulars of the band structures are as follows: 

In the NaC£ structure the gap is always indirect, between the top of the 

chalcogen valence p band at r and the bottom of the Ba 5d band at X. The 

shape of this d-band is similar to that in a fee noble or transition metal, 

* 

as would be expected from the fact that the Ba atoms constitute a fee sub- 
lattice. The Ba 6s level, doubly occupied in the atom, is empty in the 
solid because the two 6s electrons are transferred to the chalcogen valence 
p band. The results for the CsC£ structure are similar, except that at 
zero pressure the top of the valence band is at M and the bottom of the 
conduction band at r. 

Little quantitative significance can be attached to the calculated 
band gaps in the NaCi structure, which are too sensitive to the fraction 
fBa of the unit cell filled by the Ba sphere. The value of fjja for 
each volume and crystal structure is chosen by approximately minimizing 21 
the calculated total energy with respect to fga* This is a more appro- 
priate variational principle than, for example, minimizing the band gap, 
since density functional theory must obtain the correct total energy but 
not necessarily the correct band gap. As seen in Table I, a change in 
fg a of 0.U5 in either direction from the starred value (the value closest 
to minimizing the total energy), which corresponds to a change of only 6% 
in the Ba sphere radius, can change the gap for BaTe in the NaCz structure 
by as much as 0.2 eV. The difference between results obtained using an spd 
basis set and the spdf basis set described above, in this and subsequent 
tables, illustrates the effects of the Ba 4f levels. These should be 
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signif leant because of the large atomic number of Ba. As seen In the 
results, the 4f effects are small but not negligible. Unless otherwise 
noted, all stated results are obtained using the spdf basis'. 

In contrast to the NaC* case, the corresponding changes in the CsCa 
structure gaps are practically negligible. The higher coordination number 
in the CsCA structure provides better justification for sphericalizlng the 
Wigner-Seitz cell. We may therefore consider the results obtained for the 
CsCa structure to be a fairly accurate test of local density theory rather 
than of the particular approximation used to solve the effective 
one-electron Schrodinger equation. To obtain a similar test for the NaCt 
structure, calculations must be performed using a larger number of ASW 
spheres or more sophisticated (and computationally costly) band structure 
techniques, such as the linear-augmented- plane-wave method 17 . 

Despite the sensitivity of the NaCA structure band gaps to f|j a , a 
reasonable upper bound to the band gap produced by solving the one-electron 
equations exactly might be obtained by maximizing the gap with respect to 
fBa* (although we emphasize that there is no justification for believing 
that this procedure gives the correct band gap). The resulting gaps 

are compared with the experimental values 22 in Table 2. The 

g.max r r g 

theoretical gaps are low by 50-60%, a discrepancy larger than the 30-50% 
error usually associated with local density band gaps. However, the exper- 
imental trend of increasing band gap with decreasing atomic number is re- 
produced by the calculations. Furthermore, the absolute magnitude of the 
error, roughly 2 eV for all three compounds, is not particularly large when 
compared to results for similar compounds. Recent band calculations 14 * 16 
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for Csl, a compound isoelectronlc to BaTe, have yielded band gaps too small 
by 2 1/2 eV, with a fractional error of 40X. 


t h 

The last column of the table shows the gaps Ej * at the values of fg a 

which minimize the total energy. All are much smaller than the 
experimental gaps. 

We now turn to the band overlap transitions and the NaCt-CsCi 
structural transitions. The volume dependence of the band gaps In the CsCi 
structure is shown in Fig. 1. All of the gaps drop monotonlcally and 
transform from direct to indirect under compression, except for BaTe which 
remains indirect. The chemical trend in the gaps is unaffected by the 

compression. The calculated values E^ h ^CsCl^ the band 9 a P s 


CsCt structure at the volumes 2 * 3 ’ 4 v £sCi C - ' m,Tiediate ^ after the structural 

transition are given in Table 3. All of the gaps arepositive, indicating 
that the band overlap transition occurs after the structural transition. 
This is consistent with recent band gap measurements 4 on BaTe and BaSe and 
previous visual inspections 2 of BaS in the CsC l structure. The volumes 



at which the band overlap transitions occur are given in the third 


column of the table. As expected, the amount of compression required to 
induce the transition increases with the zero pressure band gap (cf. Table 


2 ). 


Although experimental measurements of 



have not yet been made, it 


is useful to compare our results with the volumes V d ^ obtained from a 
simple dielectric theory 1 originally due to Herzfeld. As seen in the 
Table, the band theoretic results. are consistently higher than the 
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dielectric theory estimates, by 10-15%, although both estimates show the 
same chemical trend. Band calculations 14 * 1S * 16 for Csl, Xe and Ar have 
predicted band overlap transition volumes higher than the dielectric theory 
estimate by roughly the same amount. 

At this stage our starting belief that the metallization volume (and 
hence pressure) would be well predicted is somewhat battered. As shown in 

Table 4, the calculated equilibrium volumes In The NaCi structure are 


16-17% lower than the experimental values 23 The results for the 

CsCi structure are expected to be more accurate than the NaCt results, 
because as mentioned earlier the larger coordination number provides more 
justification for sphericalizing the Wigner-Seitz cell. However, as seen 


in Table 4, the calculated equilibrium volume V 


th 

CsC l 


in the CsCt structure 


is still roughly 10% lower than an experimental value obtained 4 by 
backwards extrapolation of the equation of state measured beyond the 

structural transition. Since the difference between and In 


t h 

BaTe is less than the error in Vj^ , The calculated pressure must be off 

by more than a factor of two in the harmonic approximation, and probably 
much more when anharmonic effects are included. Nevertheless, we hfve 
displayed the calculated pressures in the caption to Fig. 1. To obtain 
accurate estimates of the transition pressures in the Ba chalcogenides, it 
will be necessary to devise one-electron potential which give more accurate 
equilibrium volumes than the local density approximation. 
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Table Captions 

Table 1. Sensitivity of calculated band gap in BaTe to the fraction fg # 
of unit cell volume occupied by Ba sphere, for fixed volume 
3 

(578 a^/atom pair) in both NaCt and CsC* crystal structures. 
Starred volumes minimize the total energy. 


Table 2. Comparison of calculated and experimental band gaps in NaCt 

th 

structure. E obtained using the sphere radii which minimize 
9 o* 


the total energy, 


Table 3. Quantities relevant to structural and band overlap transitions. 


V-Pt i S the measured equilibrium volume in the NaCt structure. 


The Herzfeld criterion 1 predicts the metall .zatlon volume V d,e ^ 
at which the dielectric constant is infinite in terms of the 


expt 

zero-pressure dielectric constant c and experimental volume V N q 

by (c-l)/(c+2) ■ Other quantities are defined in 

the text. 


Table 4. Calculated and measured equilibrium volumes In the NaCi and CsCt 
structures. v *aCi 9 1ven 1n a o^ atom pa1r * 
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Table 1 



• 


NaCi t 

- 

CsCi 

f Ba 


E g (eV) 


E g (eV) 


spdf 

spd 

spdf 

spd 

0.25 

0.75 

0.64 



0.30 

0.92 

0.82* 



0.35 

1.12 

1.01 

0.73 

0.64 

0.41 

1.30 

1.18 

0.72 

0.64* 

0.48 

1.40 

1.28 

0.71 

0.65 


i 

* 

f 
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Table 2 


^th.max 

9 


(eV) 



E*\(eV) 
y * 


BaS 

2.1 

3.9 

1.6 

BaSe 

_.8 

3.6 

1.3 


3.4 


0.9 


BaTe 


1.4 
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Table 3 



th 

E 9 

, st rue 
< v CsC£ ) 

yStruc 

V CsCt 

yexpt 

NaCi 

V B0 

V CsCt 

yexpt 
NaC je. 

ydiel 

yexpt 

NaCi 


spdf 

spd 


spdf 

spd 


BaS 

1.0 

0.8 

0.77a 

0.60 

0.62 

0.55 

BaSe 

0.8 

0.6 

0.78b 

0.63 

0.66 

0.58 

BaTe 

0.2 

0.0 

0.78b 

0.74 

0.77 

0.63 


a Ref. 2 b R e f. 4 
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Table 4 
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expt 

V 

yth 

NaU 

yexpt 
CsC i 

yth 

CsC* 



NaU 

yexpt 

NaU 

yexpt 

NaU 

yexpt 

NaU 




spdf 

spd 


spdf 

spd 

BaS 

440 

0.83 

0.84 

-- 

0.76 

0.78 

BaSe 

483 

0.84 

0.85 

— «. 

0.79 

0.80 


BaTe 


578 


0.84 


0.84 


0.91 


0.80 


0.81 
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F1 gure 
Fig. 1 


Caption 

. Volume dependence of band gaps of Ba chalcooenides in the CsCji 
structure. is the measured equilibrium volume in the NaCi 


structure (Ref. 23). Filled circles denote calculated points. 
Pressures at band-overlap transitions are 320 Kbar for BaS, 170 
Kbar for BaSe, and 40 Kbar for BaTe. 





